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One of the important findings at RHIC is that high transverse momentum {pr) hadron pro- 
duction in central Au+Au collisions is suppressed compared to (properly scaled) p+p collisions 
[|l|, ^. This suppression can be attributed to energy loss of high-/7r partons that traverse the hot and 
dense medium formed in these collisions. An important goal of the study of heavy ion collisions 
is therefore to determine the medium density from the measured suppression by comparing the 
measurements to sufficiently detailed calculations of parton energy loss and hadronisation in the 
medium. 

Four different theoretical descriptions of radiative energy loss in hot and dense QCD matter 
are generally used. Two of the descriptions: the opacity expansion and multiple soft scattering ap- 
proach, have been shown to be different approximations of the same basic path-length formalism 
which was formulated by Baier, Dokshitzer, MuUer, Peigne and Schiff (BDMPS) ^ and indepen- 
dently by Zakharov [Q]. The higher-twist formalism was formulated by Wang and Guo ^ and 
starts from higher-twist terms which are closely related to higher twist in Deeply Inelastic Scatter- 
ing calculations. It is also possible to calculate energy loss in Hard Thermal Loop field theory, as 
explored by Arnold, Moore and Yaffe (AMY) [^. 

For a realistic calculation of final state hadron suppression, the theoretical description of the 
fundamental energy loss process needs to be supplemented with parton production and fragmen- 
tation and a realistic collision geometry. Existing calculations have used various simplified ge- 
ometries, including fixed (mean) path length, (uniform) hard spheres and Woods-Saxon overlap 
geometries. Some recent work has been done with realistic hydrodynamic media 

For a systematic evaluation of the medium density and the uncertainties on the density, we 
would like to compare the various theoretical frameworks using the same geometry and then vary 
the assumptions about the geometry in an independent way. Unfortunately, the existing literature 
does not provide such a systematic comparison. 

For example, a long-standing issue is the comparison of the opacity expansion and the multiple- 
soft gluon emission approximation. The basic energy loss kernels for both formalisms have been 
compared in some detail and it was found that both formalisms could be made to roughly agree 
JTO[], but no connection was made to the medium density or temperature. More complete calcu- 
lations of hadron suppression using the two formalisms have been performed by two independent 
groups, which makes it difficult to compare the results. For the multiple-soft scattering approach, 
the PQM model implements a Woods-Saxon geometry where the medium density is assumed to 



scale with the local binary collision density Ql IQ. It was found that the average transport coefficient 
q that a parton sees is 13.2 GeV^/fm [1I2|], corresponding to an estimated average medium temper- 
ature (see section [O] below) of 1 GeV. A similar calculation for the opacity expansion (GLV) was 



performed by Wicks, Horowitz, Djordjevic and Gyulassy []13|], assuming the medium density to 
be proportional to the local participant density. In this calculation, it is found that a medium with 
about 1400 soft gluons in the transverse plane at mid-rapidity is needed to reproduce the measured 



nuclear modification factor /?aa[121. This value corresponds to a temperature of approximately 
370 MeV. From this, it seems that medium densities extracted using the two different formalisms 
are vastly different (a factor 2 in T corresponds to a factor 8 in density). In the following, we will 
explore some of the possible origins of these differences in detail. 
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1. The TECHQM brick problem 

A general calculation of energy loss effects on high-7?r spectra in heavy ion collisions involves 
not only the energy loss kernel which is calculated in (perturbative) QCD, but also the input parton 
spectrum, parton fragmentation, and the geometry of the reaction zone. As a result, it can be 
difficult to judge whether differences in e.g. the obtained medium density/effective temperature are 
generated by the energy loss formalism itself or some of the additional assumptions. 

The 'Theory Experiment Collaboration on Hot Quark Matter' (TECHQM) has proposed a set 
of benchmark problems to compare different energy loss formalisms, also known as the 'brick 



problem' ||14|]. The TECHQM Brick Problem asks to calculate quark energy loss for a fixed-length 
homogeneous medium (the brick). Calculations are performed for two typical path lengths: 2 and 
5 fm and 2 partons energies 10 and 100 GeV. The following discussion will be largely based on 
this problem. 

1.1 Defining a common scale: T 

In the following, we will focus on comparing the first order opacity expansion (GLV, single- 
hard scattering) [ [l5| ] approach and the multiple-soft (ASW) [ [Io| ] approach. Energy loss in any 
scheme is expected to be governed by the medium density p and the path length L. In the two 
schemes considered here, three variables are generally discussed: the average momentum kick per 
radiated gluon {kr), the mean free path A and the path length L. For the multiple soft scattering 
approach, the transport coefficient q = {kj) /X is also often used. 

In the following, the temperature T will be used as a common scale. The density of gluons in 
a QGP (or pure gluon gas) is 



For simplicity, we approximate the medium as a pure gluon gas. The gluon-gluon transport cross 



section Oog (to logarithmic accuracy) is [16] 



971 

2inL ' 



Ogg = ^, (1-2) 



where a, is the strong coupling constant and mj) is the Debye screening mass (mj) = \fA%aJ^ for 
a pure gluon gas). The mean free path is 

We further take (A:^) = m^. This ties all the parameters in the calculations to T . 

2. Comparing energy loss - importance of distribution width 

Figure |I], left panel, compares the probability (density) distribution for energy loss P{AE) 



for the WHDG opacity expansion (black curve) [13] and ASW multiple soft scattering approach 



(red curve) [|l^] for a 10 GeV quark and L = 2 fm. For this comparison, the medium densities 
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Figure 1: Left panel: energy loss distribution for a quark with E = 10 GeV in a homogeneous medium 
of L = 2 fm and (AE/E) ~ 0.2, from the multiple-soft scattering ASW approach and the WHDG opacity 
expansion. The points on the left and right axis give the discrete probabilities to lose no energy and to lose 
the entire energy. Right panel: Nuclear modification factor for pions from quark fragmentation, using the 
energy loss distributions in the left panel and KKP fragmentation. 

have been set to obtain an average fractional energy loss of 0.2 for both models (T = 500 MeV and 
q = 3 GeV^/fm). The two different approaches give significantly different energy loss distributions, 
starting from the probability to lose no energy at all P{0), which is 0.5 for the opacity expansion 
and 0.6 for the multiple soft scattering approximation. The continuous part the distribution shows 
correspondingly larger probabilities for WHDG than ASW. In addition, the opacity expansion curve 
is more strongly peaked towards small energy losses. 

A mean fractional energy loss of 0.2 was chosen for Fig. |l|, because for a power law spectrum 
with « = 8, this is expected to generate Raa ~ (1 — 0.2)"^' = 0.21 close to the measured value of 
0.2-0.25 at RHIC [[I], ^, |T2|]. The right panel of Fig. |l| shows Raa calculated for quark fragments, 
using a realistic (LO pQCD) parton spectrum and KKP fragmentation |[T7[]. The observed values 
of 7?^ is much larger than the naively expected value of 0.21, due to the broad distribution P{AE). 
In addition, it is seen that Raa is smaller for the opacity expansion formalism than in the multiple- 
soft scattering approximation, as one would expect from the smaller value of P{0) (Raa has to be 
greater than or equal to P{0)). 

A better estimate of Raa is the weighted average: 

Rn = f' de{l-eYP{e), (2.1) 

JO 

where e = AE/E. The calculated values of R^ are shown in the figure and are in approximate 
agreement with 7?^. 

3. Single-hard and multiple-soft scattering 

Figure ^ shows a comparison of energy loss in the opacity expansion and the multiple soft 
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T (GeV) T (GeV) 



Figure 2: Mean fractional energy loss (left panels) and (right panel) as a function of medium temperature 
for a quark passing through a homogeneous medium of length L = 2 fm (upper panels) and L = 5 fm (lower 
panels). The red curves shows the result in the multiple-soft scattering approach (ASW)[p^ and the blue 
curves are for the opacity expansion (WHDG)[pj|]. 

scattering formalism as a function of the temperature of the medium. The upper panels show 
results for L = 2 fm and the lower panels for L = 5 fm. The values of q for the ASW calculation 



and /I, A for the WHDG calculation were calculated using the equations in Section |l.l| . It can be 
clearly seen in the figure that the WHDG formalism typically leads to a larger mean energy loss 
or smaller Raa for a given temperature. This might be due to the larger impact of the high-energy 
tail in the energy spectrum of radiated gluons in the opacity expansion than in the multiple soft 



scattering approximation []10|]. The differences are more pronounced at low temperatures and short 
path lengths than at high temperatures and long path lengths. For L = 2 fm, the fractional energy 
loss reaches 0.2 at T = 500 MeV and T = 600 MeV for WHDG and ASW respectively. The 
corresponding values of are close to 0.6. With a medium length L = 5, the energy loss is much 
larger, reaching ~ 0.2 at T = 375 MeV for WHDG and T = 530 MeV (<? = 1.9 GeV^/fm) for 
ASW. Clearly, large path lengths are needed to reach the value Raa ~ 0.2 that is observed at RHIC. 
The corresponding mean energy loss is also large: £ ~ 0.5. 

It is interesting to note that using the brick geometry, the difference between the temperature 
needed to reach Raa = 0.2 in the GLV and BDMPS formalism is about a factor 1.5, while some 
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Figure 3: Distributions of effective path length L^fj and transport coefficient qf,ff in the case where the 
local medium density is proportional to the density of hard scatterings p^aii (solid line) and the density of 
participant nucleons Ppan (dashed hne) for 0-10% central Au+Au collisions. The proportionality constants 
for the medium density scaling are given in the text. 



earlier studies indicated factors of 2-2.5 (see introduction). 



4. Note on geometry 



The space-time geometry of the collision is an important ingredient of realistic calculations of 
hadron suppression due to parton energy loss. In the above, a fixed path length and homogeneous 
medium were used to set up a simple benchmark comparison. Most existing calculations in the 
literature use somewhat more realistic, parametrised geometries which start from hard spheres 
or the Woods-Saxon density profile of cold nuclei. Only very recent calculations use realistic 
geometries from hydrodynamic calculations [Q-^. 

In particular, the PQM model [11] uses quenching weights from the multiple-soft scattering 
approximation with a density profile that is proportional to the local density of hard scatterings 
Pcoii, while the WHDG calculations [13] use the opacity expansion with a density profile that 
is proportional to the local density of participants Ppan- Figure |3| compares the distributions of 
effective path lengths L^ff and transport coefficient q, for a medium medium density scaling with 
Pcoii and Ppart- l-eff ^nd q^ff are calculated using the same averaging procedure as in the PQM 
model: 

with 

poo 

duu"q{u), (4.2) 

where q{u) is the local density and the integral runs along the path of the parton through the 
medium. The points of origin for the partons follow a p^oii distribution in both cases. 
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The density distributions as a function of the position the transverse plane s and the impact 
parameter b are calculated from the thickness functions Ta using 

Peoii{s;b) = Ta{s- \b) Tb{s + ig), (4.3) 

Ppart{s;b) = Us-\b) (^l-e-''^'^^'+l^^'^+Ts{s + \b) (^l-^.--^'»(-5^)). (4.4) 

The terms in brackets in the equation for Ppan are needed to cut off the naive Ppart °^Ta + Tb in 
cases where the density in one of the nuclei is close to zero. These terms are the Poisson probability 
to have at least one struck nucleon in the 'other' nucleus (a is the total inelastic cross section for a 
proton-proton collision, 42 mb at RHIC energies). 

There is one free scaling parameter in both distributions. For the case of Pcoii scaling, the 
parameter has been set to the same value as in the PQM model {kcou = 5 • 10^ fm/A^ \^^)- The 
value for Ppan was chosen such that the average qe/f was similar to the Pcoii case {kpan = 2.3 • 
lO^^fm-V^)- 

The qualitative difference between the Ppan geometry and the Pcoii geometry is that the par- 
ticipant density is more uniform than the collision density. This leads to a larger L^ff and a more 
sharply peaked distribution for q^ff than in the case where the density scales with Pcoii- The cor- 
responding suppression is therefore expected to be larger for a participant density geometry than 
for the collision density scaling. This may partially explain the differences between the medium 
density as determined from the WHDG model and from the PQM model. 

5. Discussion and outlook 

A full calculation of hadron suppression due to parton energy loss in heavy ion collisons 
contains a number of ingredients. Existing calculations have found a large range of different values 
for the medium density that is needed to describe the measured suppression at RHIC. 

When comparing radiative energy loss for a fixed pathlength in a homogeneous medium in the 
first order approximation of the opacity expansion and the multiple soft scattering approximation, 
the differences are found to be sizable, with a significant dependence on path length and medium 
density or temperature. For example, Raa ~ 0.2 is reached for a quark in a gluon gas of temperature 
T = 375 MeV in the opacity expansion and T = 530 MeV for the multiple soft scattering approach. 
In both calculations, the average energy loss is large (approx. 50%) and the probability distribution 
is very broad. This is an important to keep in mind for phenomenolgy, because it means that one 
cannot think in terms of 'typical' or 'average' energy loss. In addition, the formalisms to calculate 
energy loss use approximations that are applicable for small fractional energy loss. It would be 
interesting to estimate the uncertainties from these assumptions. 

Another important aspect of energy loss models is the medium density profile. Two different 
assumptions are commonly found in the literature: the local medium density is either proportional 
to the local participant density Pp^rt or to the local collision density Pcoii- The latter approach leads 
to shorter effective path lengths and larger fraction of partons that experience a small effective 
medium density. Both effects lead to a smaller energy loss for a Pcoii density profile than for a Ppan 
profile with the same central density. 
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Future work will concentrate on using the different energy loss formalisms in a more realistic 
geometry. We also aim to extend the comparison to the Higher Twist and Hard Thermal Loop 
formaUsms. 
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